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We study the dynamics of classical particles in different classes of spatially extended self-similar 
systems, consisting of (i) a self-similar Lorentz billiard channel, (ii) a self-similar graph, and (iii) a 
master equation. In all three systems the particles typically drift at constant velocity and spread 
ballistically. These transport properties are analyzed in terms of the spectral properties of the 
■ operator evolving the probability densities. For systems (i) and (ii), we explain the drift from the 

' properties of the Pollicott-Ruelle resonance spectrum and corresponding eigenvectors 

(N ' 

03 
i— s 



- 1—^ 

X 
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I. INTRODUCTION 



Billiard systems have long served as paradigm models to study the foundations of statistical mechanics in connection 
to ergodic theory 0,0. In the recent years the study of transport properties of ensembles of particles in spatially 
extended systems, like diffusion in the Lorentz gas or multi-baker map 0, 3 and heat conduction in similar systems 
0- E3- 01 - has proven to be very fruitful in establishing connections between irreversible phenomena at the macroscopic 
scales and the chaotic properties of the reversible classical dynamics at the microscopic scales 0, |t| . 

In this respect an essential tool is the Liouvillian formulation of the dynamics. In this formulation, instead of 
considering the behavior of individual trajectories, we consider the evolution of a density po{X) of initial conditions. 
1 — This density evolves in phase space according to the Liouville equation d t p t (X) + Lp t (X) = 0, where Pt(X) represents 
I . the density at time t in phase space of points X and the operator L — {H,-} with H the Hamiltonian of the system 
' and {•, •} the Poisson bracket is called the Liouvillian operator. This equation is integrated using the initial condition 
(N| . Po(X) ^(|- We write its solution in the form pt(X) — (P t po)(X) where we introduced the evolution operator P* 
known as the Perron-Frobenius operator. When we are interested in the future time evolution of the system we may 
analyze this operator in terms of the Pollicott-Ruelle resonance spectrum {sj} of the chaotic systems ^lj, which 
determines the decay rates of the system. That is, for long times, the density can be decomposed on modes which 
decay exponentially in time p(X,t) ~ e Sjt Cj?pj(X) with the Cj determined by the initial condition. Within this 
■ theoretical framework, macroscopic properties have been related to microscopic quantities. In particular it has been 
""^ I applied to billiard systems which are spatially periodic and whose extension is infinite in one or two directions y|. 
For these systems, analytical results can be obtained using the Bloch theorem. In fact, expanding the functions 
like densities p and eigenstates of the Liouville operator ipj in Fourier series it is possible to analyze the problem in 
the finite domain of a unit cell instead of the infinite domain of the extended billiard, and compute quantities like 
^ ' eigenvalues as function of the wavenumber. 

In studying the application of this formalism to situations of physical interest, an important problem concerns the 
characterization of classes of systems other than spatially periodic ones, where one can successfully apply the spectral 
theory of the evolution operator. 

Here we study a class of extended billiard models with a self-similar structure and show how the techniques described 
above can be transposed to such systems. By self-similar billiard, we mean a billiard made up of a collection of cells 
with a one-dimensional lattice structure, where the cell sizes increase exponentially with their indices. Because of this 
property, particles will move in a preferred direction and therefore have a mean drift. That is, contrary to the periodic 
case, a density of particles drifts with constant velocity and does not diffuse. It is our goal to provide a theoretical 
understanding of these properties, based upon the spectral analysis of the evolution (Perron-Frobenius) operator of 
the system. 

By understanding we mean that, although we do not obtain explicit solutions for the eigenvalues and eigenvectors 
of the evolution operator, we show that they verify two properties which are essential to produce this drift. Thus, we 
establish the connection between a transport property of macroscopic nature and the evolution operator acting on 
phase-space trajectories, for a new class of spatially extended systems. 

A comparison between different levels of description is achieved by considering successively a fully deterministic 
self-similar billiard, then introducing a mesoscopic model with stochastic collision rules, and finally a marcoscopic 
model in the form of a master equation. 
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The article is organized as follows. Section [H] describes the self-similar billiard and the evolution of a particle in it. 
We identify two properties which characterize the spectrum of the evolution operator and discuss their consequences 
on the evolution of statistical ensembles. In order to help understand these features, we introduce, in Sec. IIIII a 
class of self-similar graphs. We show this system verifies two properties similar to those of the billiard. In Sec. IIVI 
a phenomcnological approach is given, based on a Master equation, for which we obtain exact expressions for the 
drift velocity and the mean square displacement. This provides theoretical predictions for the billiard and the graph 
models which are compared to numerical computations. Conclusions are drawn in Sec. Ivl 



We will consider self-similar billiard chains such as shown in Fig. |2 which consist of an infinite collection of two- 
dimensional cells, shown in Fig.^ glued together along a horizontal axis. Each cell contains convex scatterers and is 
open so as to allow particles to flow from one cell to the next. The shapes of the cells are identical, but their sizes 
are taken to grow exponentially with their indices. The overall geometry is such that upon combined shifting and 
rescaling the whole billiard is unchanged. 



We consider a self-similar billiard chain based on the Lorentz channel y| • The reference cell is represented in Fig. ^ 
It is a region defined by the exterior of five disks, four of which are half disks, located at the corners of the cell and 
shared with the neighboring cells, and one located at the center of the cell. The dissymmetry between the left- and 
the right-hand sides depends on the scaling parameter /i (fi = 1 is the symmetric case). Given the value of /i, there are 
three other parameters, namely D, R and r, which, as shown in Fig.^ determine respectively the horizontal width of 
the reference cell, the radii of the external disks and the radius of the center disk. Of these three pamarameters, only 
the ratios R/D and r/D are actually relevant. Notice the mirror symmetry of the unit cell about its center under 
the transformation [i— > Appendix 1X1 details the restrictions imposed on the values of the parameters, chosen so 
that the self-similar billiard shares the hyperbolicity of the Lorentz channel. 

The two vertical segments of lengths A/^/JI and Ay/Ji, with A = \/3D — 2R, at the left and right boundaries will 
be referred to as the windows of the cell, because a particle that goes across them moves from one cell to one of its 
neighbors, as will be detailed below. 



FIG. 1: Left Panel : The reference cell geometry. A detailed description is presented in appendix [X] Here it is sufficient to 
note that A = \/3D — 2R and that with this choice the usual Lorentz channel is retrieved for fj, = 1. Right Panel : The arrows 
on the right pannel represent schematically the matching condition : A trajectory leaving a cell from the right is reinjected to 
the left with a change of vertical coordinate y = h —> y = ^ and the velocity changed according to v — > ^ . 

It will be convenient to take D = 1 and rescale R and r accordingly. Thus our billiard is characterized by fx, R/D 
and r/D. 



II. SELF-SIMILAR BILLIARDS 



A. 



Definition of the Model 
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FIG. 2: The self-similar billiard and a trajectory for fi > 1, 1/fi < 1 and [i — 1 respectively. Notice the symmetry of the billiard 
under \i — > and x — > —x. 

The whole chain is constructed by adding a cell to the right of the reference cell, identical in shape but with all its 
lengths multiplied by /x and another one to the left with all the lengths divided by /j,. We repeat this construction 
in such a way that in the ith cell to the right all the lengths are multiplied by [i % and, equivalently, by to the 
left. The resulting billiard chain depicted in Fig. [21 is so constructed that the mirror symmetry with respect to the 
transformation fi — > remains. 

Now we consider a particle moving inside the billiard with velocity v. Figure [2] shows such a trajectory. As the 
particle moves from one cell to a neighboring one, the length scales change by a factor /i, so that the characteristic 
time between collisions with the walls changes accordingly (the speed stays constant). Equivalently we can rescale 
the velocity by /i and keep the length scales unchanged. That is, in going from one billiard cell to the next, say from 
left to right, the following transformations are equivalent : 




In both transformations the time between collisions with the walls is shorter (/.i < 1) or longer (/.i > 1) by a factor \i. 

Therefore we can analyze the dynamics on the self-similar billiard in terms of the dynamics in a periodic billiard 
if, instead of rescaling the size of the cell, we rescale the velocity. This way, the dynamics on the infinite self-similar 
billiard chain can be reduced to the dynamics on a single cell, provided the matching conditions given in Tableware 
imposed. 





Exit to the right 


Exit to the left 


position r : 
velocity v : 


(D,h) -» (0,h/fi) 

{Vx,Vy) {Vx/H,Vy/n) 


(0,ft') {D,h'n) 

{Vx,Vy) -> (v X fJ,,Vy(J,) 



TABLE I: Matching conditions for trajectories crossing from one side of the cell to the other, as depicted on the right pannel 
of Fig. Q The notations h and h! are detailed in the text. 

In Tabled h is the y coordinate of the trajectory escaping through the right window, — < h < *//lA/2 and 
h', that of the trajectory escaping through the left window, — A/(2^//I) < h' < A/(2 % //7). For the collisions with the 
walls the dynamics is determined by the Birkhoff map 8j (s, Vt) n (s, v t ) n +i of our modified Lorentz channel. Here 
the variable s represents the arc-length along the unit cell boundary and Vt the projection of the normalized velocity 



4 



(\v\ — 1) to the vector tangent to the boundary. This map, together with Tabled provides the map for the evolution 
of a particle in the self-similar billiard, using only one cell. The matching conditions given in tableware analogous to 
periodic boundary conditions for a periodic billiard without the self-similar structure of our billiard. 

The construction of the self-similar billiard we have presented is general and can be used to construct other self- 
similar billiards by changing the choice of the unit cell. The crucial point is that the cells are scaled uniformly by the 
factor fj, at every step of the hierarchy. 

To our knowledge, we are the first to consider a billiard with such a geometry. 



B. Poincare Map 



As we have shown, the evolution on a self-similar billiard can be considered on a single cell, provided we change 
the speed of the particle at every time it crosses the windows. The mapping from one point of the boundary to the 
next can be described by the variables s and v t . Let £ denote this pair of variables. The map, 

£n+l = <K£n), (2) 

which determines the sequence of points visited on the boundary and the corresponding projection of the normalized 
velocity to the tangent vector, is called the "Poincare map". The area where £ lives defines the Poincare surface of 
section V '. 

This Poincare map misses the information on the speed, which we can restore as follows. To consider the change 
of speed we need to keep track of the cell where the particle is located after the n iterations. Let us introduce a new 
variable /„, which takes integer values and labels the cell where the particle is at the nth iteration of the map. We 
defined the jump function a(£), such that 

In+l = In + a(£ n )> (3) 

where a(£ n ) = 1 if 4>(Cn) has the spatial coordinate s„+i on the right window, and a(£ n ) = ~ 1 if </>(£n) has the spatial 
coordinate s n +l on the left window. Otherwise a(£„) = 0. 

Using the variable / we can determine the actual speed from Eq. Q. Equivalently, we can say that the time it 
takes a trajectory at a (phase space) point £ on the boundary to intersect again with the boundary of the billiard 
depends on / as 

T(£,/) = ^// (4) 
v 

where L(£) is the length of the trajectory between intersections at £ and <f>(Q of the trajectory with the boundary of 
the unit cell and v is the speed of the particle. 

Now we have a complete description of the dynamics in the self-similar billiard. Every point of phase space is 
characterized by X = [£,r, I] with < r < T(£, /) a new variable that restores the position between collisions. The 
complete flow <J>*(X) can be specified using these coordinates. This decomposition of the flow in terms of a Poincare 
map and the variable r along the trajectory issued from the Poincare section is called a suspension flow. It can be 
implemented in billiards and other type of systems |T^ . 



C. Evolution of Statistical Ensembles 



Consider now an arbitrary distribution of initial conditions. The evolution of this statistical ensemble is determined 
by the Perron-Frobenius operator 



which requires knowledge of the backward dynamics given by 



(5) 



*-*(£, T, J) 



(£,T-t,I), if 0<t<T, 

{4>- 1 C, r - t + T(d>-^, I - a{r !£)), I - o^-'O), 
if r<t<T + T(0- 1 f,/-o(0- 1 O), 



(6) 
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At t = r we cross the section 7? and we have to identify [£,0,1] = T^" 1 ^ / - a^^Q), I - a^ 1 ^)]. The 

interpretation is the following : for a particle in the cell J, its position and velocity are completely specified by £, 
which provides both the direction of the velocity and the last point of intersection with the perimeter of the cell. 
Moving a distance r along the line going from this point of intersection in the direction of the velocity (both specified 
by £) we retrieve the exact position of the particle. 

Now running the time backwards at t = r we are just at the point of intersection with the billiard. This point is also 
the end of the segment issued from </> _1 (£), the previous intersection with the perimeter and at r = T(0 _1 £, I— a(0 _1 £)) 
which is the upper limit for the possible values of r that start at <^ _1 (£), with the corresponding direction. This 
intersection is not necessarily in the same cell I. The dependence on I — a(0 -1 £) indicate that the point <^ -1 (£) is in 
the cell I — a((/)~ 1 ^). This identification is made at every point of intersection. 

In general, we can write Eq. © under the compact form, 



*-%t,I] 



t - 1 + t\ 4>-% i - J2 . 1 - E a fo _i o 



i=l 



3=1 



(.=1 



(7) 



if o < r - 1 + ti r i t a (^ 3 o ) < T [ <t>~ n z> 1 - E a (^) 

i=l \ 3=1 J V 3 = 1 

This construction introduces a small generalization of the treatment of periodic billiards described in [l3| . More 
details can be found in that article. 

This suspension flow formalism of the billiard dynamics has two main advantages : (i) it provides a method to 
simulate the dynamics of particles in extended billiards and, (ii) it can be used (at least formally) to determine the 
spectrum of the evolution operator as shown in 0,^3- The generalization is direct, the only difference being the 
dependence of T on I, the cell index. 

Next we present the results relevant to this generalization and deduce two properties of the spectrum of the Perron- 
Frobenius operator that follow from the self-similar structure. 



D. Two Properties of the Spectrum of the self-similar Billiard 



Following [T^, by a Laplace transform of the Perron- Frobenius operator Eq. JSJ and using Eq. Q), we obtain the 
following equation for the Pollicott-Ruelle resonance spectrum {s} and the associated eigenstates {b s }, 



with R s defined by 



R a b.[Z,I\=b.[Z,I\ i 



RJ S [£, I] = exphsT^-^, I - aOT'O)]/^" 1 ^ I - 



(8) 



(9) 



In this expression, the operations must be understood in the following order : First take £ — > everywhere and 

then I —> I — a(</> -1 £). The operator R s can be considered as a reduced Perron-Frobenius operator which evolves 
densities between successive crossings of the surface of section. A formal eigenstate can be obtained by successive 
applications of this operator to the identity. That is, 



b s [£,I] 



lim R\ 

oc 

ri exp 



3 = 1 



-sT 



J|exp 

3=1 



- j '£,j-5>(r*o 

\ i=l 



which satisfies 



R s b s 



(10) 



(11) 



Here we used Eq. (0J to obtain the last equality in Eq. (|10|l . Similarly, an eigenstate b s of the adjoint operator can 
be obtained. 
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Now, we assume that a value of s, i. e. a Pollicott-Ruelle resonance, is known and thus the formal expression of b s , 
Eq. H1U|) . defines the corresponding eigenstate. We want to show that, given the pair s and (normalized) that 

satisfies 

(R,b.)[Z,I\=b.[&I\, (12) 

there corresponds a pair s/x, b s ^ such that 

(i2,A/.)K»^ = MC^» ( 13 ) 

i. e. sfi is also a resonance. 

In order to show this, we make the following ansatz : b sfjl [£, I] — b s [^,I + 1]. First notice b a ^ exists and it is 
normalized by construction. We need to check that Eq. (| 1 31) is satisfied : 

(R afl b sfi )[£,I] = cxph^Ttrti-a^O^Mrt^-a^O], 

= exp[- SM T(0- 1 e, I - air^Wsir 1 ^ I + l- ai^O], 

= exp[- S T( ( /)- 1 e, / + 1 - a^O^bsi^t 1+1- ofa" 1 *)], 

= (R s b s )[t,I + l], 

= b t [£,I+l], 

= b s ^,I) . (14) 

The first equality follows from the definition of R s , the second from the ansatz for b sfl , the third from the scaling 
property of T, the fourth using again the definition of R, the fifth from Eq. (|12f> and the sixth again from the definition 
of 6 sp . Hence s/x is a resonance and 6 S/J [£, J] = b s [£,I + 1] is the corresponding eigenstate. 

In the appendix 151 we give a proof without this assumption. 

These two properties are essential to what follows : 

Property 1 (resonances) If s is a Pollicott-Ruelle resonance of a self-similar billiard, so is s/x. 

Property 2 (eigenstates) If the eigenstate associated to s is 6 S (£,/), then the eigenstate b sfJl (^,I) associated to s/x, 
is equal to b s (£, I + 1). 

We end this section with a few remarks. These properties follow essentially from the exponential dependence of 
T(£, J) on /i and thus do not require any assumption on the geometry of the unit cell apart from those needed to have 
a spectral decomposition of the evolution operator [l^. Property says that the mode b sfl has a lifetime which is a 
factor 1/fi shorter if fi > 1 (resp. longer if /x < 1) than the lifetime of the mode 6 M , and property [21 says that, the 
mode b st j, is equal to the mode 6 M , but shifted to the left by one cell of the self-similar chain. 



E. Transport Properties of self-similar Billiard: Numerical Results 

We claim the two properties derived in the previous section govern the macroscopic behavior of a density of particles 
in the system. Indeed it will be shown in Sec. lIII Cl that these properties induce a constant drift of the particles towards 
the direction of growing cells. More information on the spectrum, such as shape of the eigenstates and values of the 
resonances s are actually needed to compute an explicit formula for the velocity, but with these two properties we 
can already understand the main behavior, i. e. the drift. Another important aspect of the macroscopic evolution is 
the spreading of particle densities around the drift. 

In order to illustrate this macroscopic behavior, let us consider an ensemble of initial conditions located in a given 
cell Iq, whose velocities are distributed at random angles, but with the same magnitude vq. We follow the evolution 
of this density and compute the average position along the horizontal axis, (X), of this ensemble as a function of 
time. In these simulations we fix two parameters, namely r/D = 0.395 and R/D = 0.480 an we vary jj, in the interval 
0.653 < fj, < 1.532 determined by the conditions of appendix lAl 

For fj, = 1 the chain is periodic and we know that the mean value stays constant (no drift) and the spreading of 
the density is that of a diffusive process. Thus the density for long times is distributed according to a Gaussian, with 
a variance growing as yi. This is confirmed by the numerical result shown in Fig. [31 central curve on panel (a) and 
bottom curve on panel (b). 

For fi ^ 1 we observe that the packet moves towards the region of growing size cells (i. e. to the left for /i < 1 
and right for fj, > 1). The velocity of this motion is obtained by computing the average position of the particles at 
different times, which is plotted in Fig. [3 Clearly there is a well defined constant speed for each value of /i. This 
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FIG. 3: (a) Average horizontal position {X) t in the billiard, as a function of time for different values of fi. We see that for each 
value of (i there is a well defined slope (average velocity) and as jj, increases so does the slope, (b) Plot of the mean square 
displacement {AX 2 )t, for the same values of fi. We observe a clear quadratic behavior ~ t 2 (upper solid line) for all values 
of /i, except for /i = 1 (lower curve) where the mean square displacement behaves linearly with time. For panel (a) the data 
corresponds to fj, = 1/1.53 (V), 1/1.5 (-), 1/1.4 (x), 1/1.3 (o), 1/1.2 (□), 1/1.1 (A), 1 (o), 1.1 (A), 1.2 (□), 1.3 (o), 1.4 (x), 
1.5 (— ), 1.53 (V). Note that these values are chosen such that we have pairs fi and 1/fJ.- The symmetry v — > —v when /i — > 1/fi 
is seen in figure (a). In figure (b) the data for fi and 1/fi almost superpose, thus we have plotted only values for /i > 1. We 
recall the parameter values are R/D — 0.48 and r/D — 0.395. 




' 0.6 0.8 1.0 1.2 1.4 1.6 

FIG. 4: Dots: Plot of the velocity v as a function of /i. The solid line is the function c(fi — 1/ /t)(y/Jl + l/\/7i), with c = 0.05 a 
constant we use to fit the data close to fj, = 1 (the reason to plot this function will become clear in Sec. II vl . The dashed line 
is a straight line with slope 0.21. Note that for values of \i away from p = 1 the dots are not well approximated by the line. 

speed is plotted as a function of /i in Fig. |3J The data is plotted together with the speed, Eq. I|35|) . obtained from a 
macroscopic model based on a master equation to be described in Sec. II VI and a linear approximation for small values 
of /i — 1. 

We also computed the mean square displacement AX 2 = (X 2 ) — (X) 2 as a function of time. The numerical 
observation, Fig. [2t> is that AX 2 ~ t 2 , i. e absence of diffusion, as expected from the macroscopic description of the 
system in Sec. IIVI 

III. A SELF-SIMILAR GRAPH 

In order to simplify matters a bit, we introduce a stochastic system in the form of a graph, with properties similar 
to the self-similar billiard described in the previous section. In a sense, this self-similar graph is an intermediate 
description between the billiard and a macroscopic description based on a Master equation, to be discussed in Sec. IIVI 

Graphs are geometrical objects made of vertices connected by bonds of given lengths, on which particles are allowed 
to move according to specific rules. The particles move freely on the bonds and get randomly scattered at the vertices 
to any of the connected bonds, according to prescribed transition probabilities. The classical dynamics of graphs share 
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many properties of one-dimensional chaotic maps, in fact every graph can be mapped to an expanding one-dimensional 
map |15| . 

The model we consider is self-similar in the sense that the lengths of the bonds grow exponentially with their 
indices. Here the chaotic dynamics of the billiard is replaced by a one-dimensional free flight and a probabilistic law 
that determines the transition probabilities between cells. 



A. Definition of the Model 



The graph is defined as follows. It consists of a linear chain of bonds, connected to one another from end to end. 
Transition probabilities to the left, Tl, and to the right, Tr, are taken to scale according to Tr/Tl = /i, with the 
parameter /i playing a role similar to the billiard's. This choice is motivated by the dynamics of the billiard. Indeed, 
provided we can assume the escape probabilities are small, we expect trajectories to distribute uniformly within the 
cells before exiting, so that the escape probabilities to the left and right are within a ratio [i of one another. Thus let 
Tl = p/ y/Ji the transition probability to the left, and Tr = Py/JI to the right. Like in the billiard, we have a reference 
cell (with index zero) with a bond whose length we take to be Lq. The first cell to the right has length L\ = Lofi and 
the ith cell to the right has length Li — L fi l . Likewise the first cell to the left has length L_i = Lq/j.^ 1 and the ith 
cell to the left L_j = Lo/-*~ 4 - The transition probabilities for all the cells are the same as for the zeroth cell and we 
take for left vertex transition Tl and reflection Rl = 1 — Tl, and for the right vertex Tr and reflection Rr = 1 — Tr. 



Tr > < T L 

<-T- T L T R r > 




FIG. 5: Schematic representation of a self-similar graph with /i = 2. 



B. Two Properties of the Spectrum of the self-similar Graph 



We identify two properties of the spectrum of the Perron-Frobenius operator of the self-similar graph, analo gou s 
to Properties 11121 for the billiard. The Perron-Frobenius operator and its spectral decomposition was studied in [lfij . 
Here we use the tools developed in that article. In particular, as explained in [15j, the Pollicott-Ruelle spectrum {s»} 
for graphs is obtained by computing the zeros of the secular equation 



det[I-Q(s)] = 0, 



(15) 



where, in the case that concerns us, the matrix Q of the infinite system has elements Q a b — Pabe sLb , given by the 
probability of going from the bond b to the bond a times the exponential function involving the length of bond b. 
Here we must take the bonds to be directed, thus b and b represent the same bond, but with opposite directions. In 
our case, the graph has a linear structure and only neighboring bonds are coupled. The matrix Q is ordered in such a 
way that the column a and row b goes in order of increasing index, alternating directions, {. . . , —1, — 1, 0, 0, 1,1,.. .}, 



Q(s) 



V 





T L e sil 









R L e sL i 








R R e sL i 








T L e sL2 


T R e s ^ 








R L e sL2 








R R e sL2 









T R e sL * 





(16) 



/ 



From Eqs. I)15 [l -I)lti [l . and from the definition of Li, we infer the following property : 

Property 3 (Analogous to property IT]) If sq is a solution of Eq. i|_?5|) then so is So£t- 

To see this, note that changing s by sfi in Eq. I|16fl amounts to shifting the elements of the matrix by two rows 
downwards and two columns to the right, because it is equivalent to changing Li — > Li + \. Since the matrix is infinite, 
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this change leaves Q unchanged. Thus, if s is a solution of the secular equation l|15(l then so is sfi. Thus, given a 
solution so, we can identify a family of solutions {silf ^ where each = s /i 4 . Given a different solution s , there is 
another family {s^} 00 ^. 

For a given family of solutions, the ratio between to roots of Eq. I|15|) is an entire power of ji. Roots from different 
families do not have this property. 

Property 4 (Analogous to property [SJ) If Xo = {■ ■ ■ 

.} is the eigenvector associated to 

So, i. e. Q(so)xo — Xo> then the eigenvector associated to sj = Sop * s Xi = {• ■ ■ i a 0: a 6> a l' a i' a 2i a 2> ■ ■ •}> whose 
elements are shifted to the left compared to xo- That is, if dj is the ith element of xo, then it is also the (i — 2)th 
element of xi ° r in general Xj+i\b] = Xj[b + 2], which implies, Xj[b] = Xo[° + 2j]. 

With the same argument we have that the left eigenvector satisfies Xj+i[fr] = Xj[° + 2]- 

Property 0| has a simple interpretation. First, note that due to the order we used for the matrix Q a shift by two 
elements on the vectors corresponds to a shift by one cell on the chain. Therefore, property 01 says that, the mode 
Xj+i is equal to the mode Xj-, but shifted to the left by one cell. And property [3] says that the mode Xj+i has a 
lifetime which is a factor l/p shorter, if p > 1, (resp. longer if p < 1) than the lifetime of the mode Xj- Thus the 
spectrum of the evolution operator of the graph has the same properties as the spectrum of the evolution operator of 
the billiard. 



C. From the Spectrum of self-similar Systems to the Drift 

As was shown in numerical simulations of particles moving in the billiard in Sec. Ill El a density of particles drifts 
towards the regions of growing cell size of the chain with a constant velocity, while at the same time it spreads 
ballistically. We will show shortly this also holds for the graph. 

A heuristic justification for the constant drift can be given thanks to the self-similar structure of the system as 
follows. Anywhere in the system, the probability is greater for particles to exit in the direction of increasing scales. 
Furthermore as the scales grow, the average time a particle spends in a cell increases accordingly. Hence we expect 
the drift to remain constant. 

Now we offer a quantitative analysis of the drift, based upon properties ^ an d 13 or, equivalently, |3] and ^ 

The spectral decomposition of the Perron- Frobenius operator P* was obtained in [TEl ITi| for graphs and in 0, ^| 
for billiards. For graphs, the time evolution of an observable A[b, Xb] has the asymptotic expression (t — > oo) 

(A) t = (A\Ppo) = ]T(A|^)e s ^|po} + • • • (17) 

3 

with 

Wi) = E / ^ S3X/v A[b',x]dx, 

ut L b' Jo 

b' — — oo 



(18) 

i 00 rh> 

&j\Po) = W E Ub'\ / e«*/*po[b',x]dx, (19) 



b' = — oo 



and 



N^^hxMxAb}- (20) 



In Eq. (|17(l the dots stand for terms whose coefficients carry subexponential time dependence which may arise from 
degeneracies of the spectrum and can be neglected in the long time limit. Similar expressions can be obtained for 
billiards with the same conclusion 0, 0] , however for the sake of simplicity we will state the results for graphs only. 

We assume the initial density to be localized on one bond, namely pt=o{b, Xb] = Sb,b $( x b — xq), i. e. all the particles 
start from x a on the bond b^. To simplify as much as possible the calculation and notation we will consider x = 
and the bond bo as the reference bond bo — 0. We are interested in the density as a function of the position and t, i. 
e. p t [b, Xb], so we define the observable A by A[b', xy\ — 8b,b'5{xb — Xy), in such a way that (A) t = pt[b, Xb], as we see 
from Eq. I|17f) . For simplicity, we also take Xb = 0, i. e. we measure the density at the beginning of the bond b. 

Using this in Eq. l(T7|. Eq. {THJ and Eq. we get 

( A )t=T' Xj[b] Wb jis[0] ' (21) 
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Now, since the spectrum is divided in families according to whether or not Sj/sk is an entire power of /z, let us consider 
the contribution of only one of these families to the sum in Eq. 121|) . 

From property^] if %j [b] is the eigenstate associated to Sj then, the eigenstate associated to Sj/fj, is Xj-i [b] = Xj [b—'A 
and corresponds to the state \j but shifted one cell to the right. It is then easy to see that Nj — [i^Nq. 

Moreover, from Property^ we have that if so is a decay rate (Pollicott-Ruelle resonance), then there is a family of 
Pollicott-Rucllc resonances associated to it given by Sj — so/i- 7 . 

Now using the expressions of the lengths = Lo[i h , we have that the contribution to Eq. (|21[1 due to this familiy 

is 

W E xAb}^e^%[0] (22) 
iV o - p 

j = -oo 

Let us split this sum in two parts, first the terms with j — —00, . . . , and second j = 1, . . . , 00. We consider for 
now the first term only, i. e. j = 0, —1, . . . , —00. 

1. For j = 0, —1, . . . , — 00, so is the largest decay rate. Therefore the component with the corresponding rate, xo, 
is the first to decay. After it has decayed, the part of the density represented by this part of the sum -which 
we refer to as the density for short- moves to the right, because, as a consequence of property all the other 
modes are shifted to the right of xo- 

2. The support of the density is shifted to the right by a distance L, corresponding to the length of the bond where 
the mode xo is centered. Let us call it b c . 

3. Then, the component of the mode X— 1 is the second to decay with a rate so/fJ- and the support of the density 
moves another bond to the right because, by property 01 the bond where the mode X— 1 is centered is the one 
at the right of b Cl i. e. it is the bond 6 c +i. This bond b c+ \ is larger by a factor /1 than b c . 

Thus if we take the distance the packet moves divided by the characteristic time (given by the inverse of the decay- 
rate), we have that, during the decay time of the first mode, the speed was Lsq and, during the decay time of the 
second mode to decay, the speed was (Lfi)so/fJ- = Lsq, equal to the first. We can obtain the same result for the decay 
of the third mode, (L/i 2 )so/M 2 = Ls Q , and so on. 

Now, it is clear that the same conclusion follows from the other terms j = 1, . . . , 00 of the sum, thus we conclude 
the packet moves at a constant speed. 

This argument shows that a constant drift of particles has to be expected in the self-similar systems as defined here. 

If we were able to obtain the exact eigenvectors we could compute numerically this sum and compare with simula- 
tions. Instead we propose a heuristic argument and assume the eigenvectors are localized around some finite region of 
the chain. A simple possibility is to assume a Gaussian shape for the left and right eigenvectors, satisfying property 

namely Xj[b] = Xj[b] = exp[ ~^~^V 2 ''^ ]■ The argument that the eigenvectors are localized can be justified with 
a perturbative calculation provided the transmission probability p <C 1, similar to |l6[ . 

For this heuristic calculation, we consider values of [i close to one, so we write fj, = 1 + e and we can compute the 



_ M) _ „_ U^l _. „2 

measures the distance from the origin to the bond b. Then we obtain 



sum in Eq. (1221) for small e. We also do the change of variable (we take Lq = 1) y = ^ * b w y + ey 2 /2 which 



(A) ~ exp[-^ - e ^ + s t(l - ey/A) - ^] (23) 

Now, we calculate the position of the maximum of Eq. (|23|) as a function of t and obtain 

2 

y m «-ey(s t/2 + 3) (24) 

We have a maximum that is linear in /i — 1, moves in the appropriate direction and is zero if fi = 1, in agreement 
with simulations to be presented next. Thus we have shown that the two properties of the spectrum are related to 
the motion of the packet. It will be desirable to have a similar understanding of the spreading of the density in terms 
of the spectrum of the Liouvillian operator. 
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D. Transport Properties of the self-similar Graph: Numerical Results 



The evolution of a density of particles in this system behaves in a way similar to the billiard chain. As with the 
billiard, a well-defined constant drift is observed, as shown in Fig. 0(a). The dependence of this drift in fj, is depicted 
in Fig. El Here we also observe that the mean square displacement AX 2 — (X 2 ) — (X) 2 as a function of time, Fig.0 
(b), behave like AX 2 ~ t 2 . This ballistic spreading of the density is also found with the macroscopic description of 
the system as given in Sec. II VI 




tV(/L tV(/L 

FIG. 6: (a) Average horizontal position {X)t in the graph vs. time, for different values of fjb. We see that for each value of 
fx there is a well defined slope (velocity) and as \i increases so does the slope, (b) Mean square displacement {AX 2 ) t , for the 
same values of fj,. We observe a clear ~ t 2 (upper line) for all values of fj, except for // = 1 (lower curve) where the behavior is 
~ t. The same values of /j as in Fig. were used. 




IV. MACROSCOPIC DESCRIPTION 



In this section we present a system which is a charicature of the self-similar systems considered so far and helps 
understand the macroscopic evolution on systems with self-similar structure. It describes the dynamics based on a 
Master equation approach F° r this model we are able to obtain analytical expressions for the (constant) drift 
velocity and the ballistic spreading of the density. 

We consider a discrete sequence of states i e Z and, associated to them, we have a conditional probability Pi(t) 
which represents the probability of being at site i at time t given some initial condition to be specified later. We will 
identify the sequence of sites with the sequence of cells of the billiard or graph. Since for the billiard and graph only 
transitions to neighboring cells are allowed, we introduce the transition rates W~ and , respectively for transitions 
to the left and to the right of the site i. The master equation rules the time evolution of the Pi(t) and for this process 
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it has the form 

d t Pi(t) = P l+ i(t)W- +1 + P l -i(t)W+ 1 - P t (t)(W 2 + + Wr), (25) 

It is a gain-loose or balance equation. To mimic the transitions that occur on the graph, we consider W~ (respectively 
W 4 ) as the ratio between the probability of going to the left, i. e. p/y/JI (respectively p^/JI for the right) over the 
characteristic time associated to bond i Lq/j, 1 /vq, i. e., 

Wr = fV l " 1/2 , (26) 
W+ = fV /2 -*. (27) 

Thus Eq. j^Sjl becomes 



_ ^0 [, -3/2 P ,.fA-L „ 3 / 2 p. , M _ („l/2 x ,,-1/2 

In matrix notation, this reads 



fti^W = T^/i -1 /x- 3/2 ^+i(t)+/i 3/2 P i -iW-(/i 1/2 + ^ 1/2 )P4*) • (28) 

ho L J 



d t P i (t)=C ij P j (t) , (29) 

where 

= ^V J >- 1/2 <Wi + M 1/2 ^,i-i - (M 1/2 + M- 1/2 )^]. (30) 
By exponentiation, since Pj(t = 0) = <5j0) 

P<(t) = cxp[£i]yJjo (31) 

Note the symmetry of £, AjQu) = C_ji(l/ [i), implies -Pj(i) — > P-i(t) upon inverting /i. 

Note that until here, the model does not contain any metric information. Only part of the self-similar structure 
is introduced through the transition rates W i . The geometrical part of the self-similar structure is implemented by 
associating a quantity Xi to each site i which is a measure of the distance from the origin to the site i. We define it 
following the geometry of the graph. If the site i — corresponds to the bond of length Lq, and the site i to the bond 
of length Li = Lq/j, 1 , Xi measures the distance from the middle of the ith bond to the middle of the 0th bond (i. e. 
the origin), viz. 

M») = ^p^ 1 -±». (32) 
2 1 — /i 

In the expressions on the RHS, we observe the symmetry Xi — > — X_i under the change [i — ► l//x. In other words 
Xi(l/fi) = -X_i(/*). 

We define the average position at time t by 

{X) t = Y J Pi{t)X u (33) 

i 

The velocity is defined by 

V(t) = d t (X) t = J2^ijPj(t) X i- (34) 

ij 

According to the symmetries above, we should have V — ► — V upon inverting /i, which we can check by direct 
calculation. We have 



= J 1 E /' '''//' ' \X, ; - Xj ) + ^(X j+1 Xj ) 

J 
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Using Eq. (|32(l . we obtain the expression 



^) = ^ £ (m-^ 1 )(m- 1/2 +M 1/2 )- (35) 

We observe that the speed has a constant value which vanishes for /i = 1 as expected. Moreover it is clearly anti- 
symmetric upon inverting [i. 

Now we evaluate the time derivative of the mean square displacement, 

d t (AX 2 ) = d t [(X 2 ) t - (X) 2 ] = d t (X 2 ) t - 2(X) t d t (X) t . (36) 
with (X 2 ) t = Ei-PiC*)^ 2 - E q- <GS3> can be written as 

d t (AX 2 ) = ^C ij P j {t)X i [X i - 2(*>t]- (37) 

y 

replacing £ given in Eq. ijSUjl in Eq. (jSTjl we get 

d t (AX 2 ) = ^ ^ M -^ J (^[ Ai -V2 (Xj _ 1 _ + x _ 2 { Xt )) + ^ 2 (X j+1 - X ){X J+l + X j - 2{X t ))], 

y 



— f /7i I 1 



(38) 



where we used Eq. (EH) and the fact that (X t ) = Vt. 



In the long time limit, i. e. t t c = ^^rr, we can neglect the constant term and therefore AX = (X ) t — (X) 



t . In other words, the spreading of the density is ballistic, similar to what we observed numerically in the billiard 
and graph. In the opposite case t<t the constant term dominates in dt(AX 2 } and the spreading is diffusive. The 
time t c — Jjj^^rj, marks the crossover from diffusive to ballistic behavior of the mean-square displacement. 

Let us compare the speed, Eq. I|35|l . obtained in this macroscopic model with the numerical results of sections III El 
and lTTTDl In the case of the billiard (see figure the agreement is poor, because the Master equation is not a good 
model for the billiard, except perhaps for /i» 1. In fact, a particle in a given cell of the billiard escapes easily in a 
few collisions and thus we do not expect that the master equation, where the cell has a uniform probability, will be 
a good quantitative model. On the other hand what is observed is that for values of/iwl the speed is proportional 
to fi — 1. In the case of the graph, see figure [7| Eq. I|35|l fits very well the data because the probability of being in a 
bond is well approximated by a uniform distribution in the limit p< 1 where many collisions occurs in average with 
the scatterers before the particle goes to another bond. Thus we expect the master equation to be a good model for 
this system. 

Finally wc show that for this macroscopic description, we also have properties similar to 2] and [21 (or [3] and and 
therefore a correspondence between the constant drift and the spectral properties of C 

The decay rates can be computed assuming a solution of the form Pi(t) = e st Qi 1 . Upon substitution of this 
expression into Eq. 1)29(1 . we find that the decay rates s must be solutions of 

det[s<% - £ij] = . (39) 

Now we proceed as with the billiard and the graph. We assume a particular solution s of Eq. (|39|l was obtained 
and analyze det[s(iSij — C{j]. Note that due to the definition of C we have Cij = fiCi+ij+i. Therefore, 

det[s^% - dj] = det[/z(s<5 i+ ij+i - C i+lfj+1 )] . (40) 

Because fj, ^ 0, we conclude that if s is a solution of Eq. I|39l) then so is s[i. 

Now, the associated eigenvector Qf is obtained from J2j[ s ^ij ~ £ij]Qj = an d QT IS obtained from J2ji s ^ij ~ 
Cij]Q s ^ — 0. Again, using the definition of C, it is easy to show that = Qf +1 . This shows we have again the 
equivalent of properties ^ and |21 (or |21 and QJ. 



1 The complete solution takes the form Pi(t) ~ Ylj e** Q% + ' ' ' i where extra terms represented by the dots are expected from Jordan 
Blocks. 
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V. CONCLUSIONS 

In this article we have studied the statistical properties of three different classes of self-similar systems and found 
that all three share very similar macroscopic properties; we observe a drift of particle densities towards the direction 
of growing scales at constant velocity, as well as a ballistic spreading of the density. 

A justification for the presence of a drift was provided, for both the billiard chains and graphs, in terms of two 
essential properties of the spectrum of the evolution operator, the first stating that for every decay-rate s, there is 
also a decay rate s/x, and, correspondingly, the second, that the eigenstate associated to sfj, is shifted one cell to the 
left of the chain with respect to the eigenstate associated to s. 

As of the Billiard, we should note that these properties hold independently of the exact shape of the unit cell; they 
are merely a consequence of the self-similar structure. The only relevant restriction regarding the shape of the unit cell 
is that the dynamics in that billiard cell must be strongly chaotic in order to have the kind of spectral decomposition 
that we assumed. 

For the self-similar graph, these properties are also rather general. Here we considered a very simple self-similar 
graph, to avoid complicated expressions, but we can take a unit cell composed of any number of bonds, say iVg, and 
by ordering the matrix Q in such a way that all the bonds of the same cell are consecutive we will obtain property 
and property 0] which will look like = Xj[b + 2iVs]. 

Let us discuss an interpretation of these two properties in term of periodic orbits. Consider the billiard. Decay 
rates and eigenstates are determined by an ensemble of periodic orbits that form a repellor . Every periodic orbit 
can be shifted by one cell to the right or to the left and generate a new periodic orbit with the only difference that the 
period is a factor of [i respectively shorter or longer. If the periodic orbits determine the support of the eigenstates, 
then it is clear that eigenstates of the same "shape" but shifted along the billiard exist and that the associated decay 
rates will also differ by a factor of /j,. The change in period (time scale) corresponds to property ^ and the fact that 
the geometry of the orbits is unchanged to property [21 

We have also provided a macroscopic description of self-similar billiards and graphs in terms of a Master Equation, 
for which we obtained an analytical expression of the velocity and ballistic spreading. Properties equivalent to ^ an d 
121 are also shared by this system. 

We note that in going from the billiard to the graph and then to the Master Equation, the stochasticity increases. 
Billiards are deterministic chaotic systems; in graphs, there is a stochastic element which acts only when the particles 
reach the vertices, while the motion within the bond is deterministic; The dynamics described by the Master Equation 
is completely stochastic with transitions possibly occuring at any time. 

We should also note that the existence of a drift is rather intuitive and simple to explain qualitatively in terms of 
the probabilities of crossing to the left or right and in terms of the time spent in the bigger or smaller cells. However, 
our aim was to provide an example where a macroscopic property like the drift of particles can be related to the 
microscopic properties of the system, i. e. the exact Liouvillian evolution operator. 

Several perspectives are open for future research, of which we now mention two. First, it will be interesting to 
characterize the spreading in terms of properties of the spectrum and, second, regarding the billiard, we would like 
to obtain further understanding of the function v(/j.). Our main observation of this work is that there is a well- 
defined speed of propagation for each value of fi. However, a determination of the speed v vs. /x is missing thus far. 
Numerical results (not all presented here) show that, as we change the parameters of the billiard, the speed displays 
a rich behavior as a function of fi. For instance, we observed functions v(n) that are monotonous with fi in some 
parameter regions, while in other regions, such as corresponding to the parameter values used in Fig. 0] we observed a 
non-monotonous function. Near /x = 1 the behavior of v(fi) is determined by the symmetry v(/j.) = —v(l/fj,). In fact, 
it follows from that symmetry, that for small values of /i — 1, we can expand u(/x) = v'(l)([x — 1 — (// — l) 2 /2) + . . . 

We expect to continue our work in these directions. In particular, we believe a better understanding of the dynamical 
system in the unit cell with boundary conditions given by Table [I] is needed, specifically regarding the properties of 
the invariant measure. Because of the open boundary condition, the invariant measure is not the Liouville measure. 
The volumes are not conserved by the matching conditions. One therefore expects fractal properties of the invariant 
measure. This will be the subject of further publications. 

APPENDIX A: FINITE HORIZON SELF-SIMILAR BILLIARD CHANNEL 

In this appendix, we extend the finite horizon condition for the usual periodic Lorentz gas to the self-similar case, 
so as to ensure hyperbolicity of the dynamics. 

In the absence of self-similar structure (fi = 1), the Lorentz channel |8( has the symmetry of a two-dimensional 
periodic Lorentz gas on a hexagonal lattice. This system is diffusive provided it verifies the so-called finite horizon 
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conditions, namely 



2R < D < AR/V3, 



(Al) 



where R denotes the common radius of the scatterering disks and D the distance between neighboring disks. The 
lower bound ensures that the disks do not overlap and the upper one that there are no free flying trajectories. 

For a mixed Lorentz gas, with alternating rows of scatterers of radii r and _R, the hexagonal structure of the lattice 
is preserved, but the corresponding channel now consists of three rows of scatterers : the upper and lower rows with 
half disks of radii R and the middle row with disks of radii r. The corresponding unit cell has a rectangular shape 
with four quarter disks of radii R at the corners and one disk of radius r at the center (where the diagonals intersect), 
with coordinates 



4 disks of radii R 



(o,±VaD/2), 

(D,±V3D/2), 
1 disk of radius r : (D/2,0). 



For this system, the finite horizon conditions become 



D > 2r, 

D > 2R, 

D < 2(r + R)/V3. 



(A2) 
(A3) 



(A4) 
(A5) 
(A6) 



Now introducing the parameter /i, we deform the cell of the Lorentz channel from a rectangle to a trapezoid with 
vertical sides scaled by 1/ ^JJi and yfJL respectively, as was shown in Fig. ^ The positions of the disks become 



4 disks of radii R : (0^/(2^), 
1 disk of radius r : (D/(l + /x), 0), 

where the disk on the horizontal axis lies at the intersection of the diagonals. 
Let us assume /i > 1 in the sequel. 



(A7) 
(A8) 



1. Non-Overlapping Disks 

The non-overlapping conditions Eqs. I|A4IA5(I transpose to three new conditions : 

1. The sum of the radii of the external disks must be less than the length of the lower/upper side of the cell, 



R 



y/JIR < 



\ 



D 2 



2. The shortest diagonal must be greater than the sum of the radii R/ \f]i and r, 

R 



(A9) 



r < 



( V3D 



D 
li + 1 



(A10) 



3. The radius of the center disk must be less than the distance from the intersection of the diagonal to the cell's 
boundary. That is, 



r < 



D 

i + n' 



(All) 
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Let p = R/D and A = r/D be two dimensionless parameters. The first of the three conditions above, Eq. I|A9(1 . can 
be transformed into a second degree polynomial in p : 



(3/4 - pV - (1/2 + 2p 2 )p + (3/4 - p 2 ) > 0. 



We must assume 



The two roots of Eq. (|A12|I are 



P< ~- 



l + 4p 2 ±2y/8p 2 -2 



3 - 4p 2 

which are real if p > 1/2. In this case, the condition (|A9|> is only satisfied provided either of 

p > p+,p < p-, 

hold. Note that p+ = 1/p-- These conditions are depicted in Fig. [S] 

Otherwise, if p < 1/2, Eq. (|A12|) has no real roots and Eq. I|A9() is always true. 

In [i 

15 
10 
5 



(A12) 



(A13) 



(A14) 



(A15) 



-5 

-10 
-15 



FIG. 8: From Eqs. ET4llAl5t . The gray region is forbidden when 1/2 < p < 
Consider the second condition, Eq. IjAlOp . which takes the form 



-(p + X^JI) 2 



(p + l) 2 +p>0. 



(A16) 



This is a sixth order polynomial in ^fp, with potentially as many roots, but is easy to solve numerically. 
With these notations, we rewrite Eq. (|A11(I as 

p<\-l. 



(A17) 



Finite Horizon 



The condition that there are no free-flying trajectories, as seen in Fig. ^\ is that the line tangeant to the row of 
upper disks intersect the disk at the center. 

Consider the reference cell and let x denote the distance from the center of the central disk to the trajectory tangent 
to the upper disks. Since these lines are perpendicular, we can write 



x = (L + y) sin /3, 



(A18) 



where y = D / (p + 1) is the distance from the boundary of the cell to the center of the disk and L = D^^p 1 
D / (p — 1). In order to compute the angle (3, we write (3 = a — 7, where a and 7 are given as follows. 
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FIG. 9: Geometry and definition of a, f3, 7 and x 



First we have that the difference between the two verticals from the base line to the upper right and upper left 
disks is v / 3/(2-D)( v / /I — 1/-JJI) and D is the horizontal distance between their centers. Therefore 



tan o = ^ ( v // - 



We note that the distance a along the wall is given by 

a = D I cos a, 
D 



= - v /3(^ + M- 1 )-2. 



We can then compute 7 from the equality 



sin 7 



2R fi-1 

D y/3f? - 2^ + 3' 



Now, since f3 = a — 7 we have 



sin/3 = sin a cos 7 — cos a sin 7, 
2p v / /7(^ - 1) 



The finite horizon condition is 



In the limit fj, — > 1 one retrieves Eq. 1A6J1 . 



A = — > 



3 - 2/i + 



2// 

D ' /i 2 - 1 



- 2 



sin/3. 



(A19) 



(A20) 



(A21) 



(A22) 



(A23) 



3. Choice of Parameter values 



We want to fix the values of p = R/D and A = r/D and find the range of values of /j, so as to verify the constraints 
specified by Eqs . (H3 lATOl IATTI |A"23|) . 

The numerical results presented in Sec. Ill El use the fixed parameters A = r/D = 0.395 and p = R/D = 0.480. 
For those values, the range of allowed values of p is 0.653 < p < 1.532. For the sake of illustration, other possible 
parameter values are shown in Figs. ITUllll with either p or A fixed. 
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FIG. 10: Values of fi vs. A, for fixed p = 0.48 (a), 0.52 (b), which are consistent with Eqs. {Kfy. tKWi . iXTTl and SKSSi . 





0.47 0.49 0.51 

FIG. 11: Values of fj, vs. p, for fixed A = 0.395 which are consistent with Eqs. <lA9t . <A"Tot . Edit and dX23t . 



APPENDIX B: ALTERNATIVE PROOF OF PROPERTIES [T] AND [U FOR A SELF-SIMILAR BILLIARD 

The dynamics in the billiard chain is Hamiltonian and therefore the Liouville measure is invariant under the time 
evolution. According to Birkhoff, the Liouville measure can be expressed as d£dr, where £ = (s,v t ) is a pair of 
variables with s representing the arc-length of the billiard (not the unit cell), and Vt the projection of the velocity to 
the tangential direction at the point s. Because the energy is conserved we can fix the speed \v\ = 1. With this choice 
— 1 < v t < 1. Now, t is a variable which measures a distance along a line drawn from the point s with the direction 
given by Vt- Obviously, r can vary from zero to the intersection of this line with the billiard. This distance will be 
denoted by T(C), *■ e. < r < T((). Thus we can consider that the average of a function /(£, r) is given by 



(/> = 



1 



fdsfdv t J^ C) dr 



no 

ds I dvt I drf{s,v t ,r). 



(Bl) 



The variable s is unbounded because the billiard is infinite. Due to the construction of our billiard, we can split the 
integral over the arc-length into the contributions over the differents cells, i. e. 



ds 



(B2) 



Due to the symmetry Lj = /j, t Lq, we can let < s < Lq. The restriction of the variable £ to this region will be called 
^. We also observe that T(£) = T(£, /) = L^)^ 1 . Thus the average can be written as 



? M 7 ds J dvt l 



drf(s,v t ,T, I), 



(B3) 
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or, in a more compact form, using J d£T(£,I) — /i / (L(^)), 



"'fi- a^w ^/*! dT/K,r,/) ' (B4) 



Now, the quantity that is of interest to us is 



1 r f T{iJ) 



lim 



We assume that s is a resonance i. e., (b s *b s ) — 1 ^| and we want to show that this implies (b lls *b tJ , s ) = 1, that is 
that s/Lt is also a resonance. 

To show this, we note that, from the formal expression of b s , we have I] = b s [£, I + 1], and therefore 

(b^b^) = lim — J 1 Yl ^ I W 8 [^ + 1]M^+1], (B6) 

which is equal to 

J+i 

(6 MS ,6 MS ) = lim — m — — ]T ^ / dtL(t)b a [t,I\b.[t,r\. (B7) 

If the limit in Eq. (|B5|1 exists (as we asume), then the limit in Eq. (|BT|) also exists and they are equal. Therefore we 
conclude that (b^b^s) = 1 and so sfj, is also a resonance. This proves property (A). 
Now, from the formal expression for the eigenstate associated to s, 

oo j 

bsltl] = Y[cxp[-sT(cb-^,I (B8) 

j=l i=l 

we can express 6 S/J [£,/], the eigenstate associated to s/i as 

oo j 

b s ^i] = n^h^(f 3 y-E fl («))]. 

°° j 

j=i i=i 
= 6aK./+l], (B9) 

where the second equality follows from the definition of T(£, /) in Eq. Q| and the last from Eq. (|B8|) . This is to say 
that the eigenstate associated to sfi is the same than the one associated to s but shifted one cell to the left, and proves 
property (B). 
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m 2/ / dZL(S)b a .{$,r)b a [S,I\. (B5) 

7 J 
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